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ABSTRACT 

EXLup is the prototype of the EXor class of young eruptive stars: objects showing repetitive 
brightenings due to increased accretion from the circumstellar disk to the star. In this paper, we 
report on medium-resolution near- infrared spectroscopy of EX Lup taken during its extreme outburst 
in 2008, as well as numerical modeling with the aim of determining the physical conditions around 
the star. We detect emission lines from atomic hydrogen, helium, and metals, as well as first overtone 
bandhead emission from carbon monoxide. Our results indicate that the emission lines are originating 
from gas located in a dust-free region within «0.2AU of the star. The profile of the CO bandhead 
indicates that the CO gas has a temperature of 2500 K, and is located in the inner edge of the disk or in 
the outer parts of funnel flows. The atomic metals are probably co-located with the CO. Some metallic 
lines are fluorescently excited, suggesting direct exposure to ultraviolet photons. The Brackett series 
indicates emission from hot (10 000 K) and optically thin gas. The hydrogen lines display a strong 
spectro-astrometric signal, suggesting that the hydrogen emission is probably not coming from an 
equatorial boundary layer; a funnel flow or disk wind origin is more likely. This picture is broadly 
consistent with the standard magnetospheric accretion model usually assumed for normally accreting 
T Tauri stars. Our results also set constraints on the eruption mechanism, supporting a model where 
material piles up around the corotation radius and episodically falls onto the star. 
Subject headings: stars: pre-main sequence — infrared: stars — stars: activity — techniques: spec- 
troscopic — stars: individual(EX Lup) 



1. INTRODUCTION 

Young eruptive stars form a special class of pre-main 
sequence objects. The class is defined by unpredictable 
brightenings of up to 5 mag at optical wavelengths, and 
it is believed that temporarily increased mass accretion 
from the circum stellar disk onto the star is responsible 
for the flare-ups (jHartmann fc KenvonllT996l ) . In spite of 
their relatively small number (around two dozen stars), 
young eruptive stars are key objects in understanding 
early stellar evolution. On the one hand, the mass in- 
flow during outbursts contributes to the build-up of the 
final stellar mass, which would otherwise be difficult to 
explain considering the observed t ypical rates of mass 
accretion (less tha n »10~ 6 M^vr" 1 . iKenvon et alJl!994t 
lEvans et al.l 120091) . On the other hand, the eruptions 
affect the density, temperature, and chemical structure 
of the disk, possibly also influencing the conditions of 
planetesimal and planet formation in the terrestrial zone 

([Abraham et al.ll2009f ). 

One of the two main groups of young eruptive stars 
is called EXors, named after the prototype EXLup. 
These objects produce smaller, 1-3 mag flare-ups every 
few years for a period of a few weeks as well as rare 5 mag 
extreme outbursts lasti ng several mo nths (see e.g. the 
light curve of VY Tau in IHerbidll977l or that of EX Lup 
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in iFfer big 2007). Although accretion is the likely source 
of power for the outbursts, fundamental aspects of the 
mass inflow process in EXors require significant clarifi- 
cation. 

For a normally accreting T Tauri star, the stellar mag- 
netic field truncates the accretion disk at a distance of 
a few stellar radii above the stellar surface. From that 
point, the material is magneticall y channeled onto th e 
stellar surface along funnel flows ([Bouvier "et~aT][2007T) . 
During outburst, the accretion rate is enhanced by a few 
orders of magnitude, and it may be the case that the 
magnetic field can no longer truncate the inner regions 
of the accretion disk and hence disk material is dumped 
onto t he star along the eq uatorial plane in a boundary 
layer (jPopham et alJU996Q . It remains an open question 
as to what path the accreting material follows in young 
eruptive stars like E Xors. 

EXLup fd=155pc. lLombardi et alj2008l ) produced its 
most e xtreme eruption in 2008 (Jones 2008; Kospa l et al.l 
l2008al) . This made it the center of attention and the 
subject of many different kinds of observations. Our 
group conducted a coordinated multi- wavelength observ- 
ing campaign at a single epoch, as well as monitoring 
studies during the outburst. 

The optical spectroscopic data of this campaign, which 
will provide insight i nto the accretion and wind proce sses, 
will be pre s ented in ISicilia-Aguilar et al.l ()in prep.f ) . In 
IGoto et all (|2011[) we analyzed the CO fundamental 
emission lines using high resolution spectra in the ther- 
mal near-infrared and found that the emission comes 
from two separate physical components: a narrow line 
component which is constant in time, and a broad line 
component which is closely related to the outburst and 
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is decaying with time. In I Abraham et al.l (|2009t ) we re- 
ported on on-going crystallization of silicate grains dur- 
ing the outburst, and our monitoring of the crystalline 
silicate features showed e vidence for a fast radial trans- 
port of silicate crystals (Juha sz et al.l [submitted]) . In 
Juhasz et al. we also modeled the whole optical-to-mm 
spectral energy distribution of EX Lup in outburst and 
found that most of the accretion luminosity is emitted as 
a single temperature blackbody of 6500 K. 

In this paper we present the medium resolution near- 
infrared spectra of our campaign, and use the observed 
emission lines to analyze the location, kinematics, and 
energetics of warm atomic and molecular gas present in 
the inner few tenths of AUs in the EX Lup system. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Near-infrared spectroscopy 

We observed EX Lup with SINFONI, an adaptive op- 
tics (AO) assisted integral field spect rograph installed 
on the UT4 telescope o f the VLT (jEisenhauer et all 
120031: iBonnet et al.l[200l ). Measurements were taken in 
service mode at three different epochs, on the nights 
24/25, 28/29, and 30/31 July 2008, as part of the 
project 281.C-5031 (PI: M. Goto). Three dimensional 
spectra were taken with the J, H and K gratings us- 
ing a Hawaii 2RG detector. The dispersion of the se- 
lected mode was 0.145 nm/pixel in the J band (1.10- 
1.45 yitm), 0.195 nm/pixel in the H band (1.45-1.85 /mi), 
and 0.245 nm/pixel in the K band (1.93-2.45 fim). The 
spectral resolution was A/AA « 2400, 4100, and 4400 
in the J, H, and K band, respectively. Spectra were ob- 
tained by rotating the instrument by 0°, 180°, 90°, and 
270°, in order to be able to better remove instrumental 
effects. Exposure time was 1 min for each rotator angle, 
giving a total exposure time of 4 min per band per epoch. 
Beside EX Lup, calibration sources of spectral type G2V 
were also observed in order to enable proper telluric cor- 
rection: HIP 93685 and HIP 79464 were observed in the 
J band on 25 July, HIP 80982 was observed in all three 
bands on Jul 29, and HIP 78652 was observed also in all 
three bands on Jul 31. 

We started the data reduction with product files pro- 
vided by the SINFONI pipeline as part of the ser- 
vice mode data delivery: product code SCDJ (full, 
coadded science product cubes) for the science tar- 
get, and PCST (full, coadded standard star cube) 
for the telluric standard. These files contain dark 
current-subtracted, flatfield-corrected, sky-subtracted, 
wavelength-calibrated, co-added data products in the 
form of 64 pixel x 64 pixel images for each wavelength, 
with a spatial scale of 12.5 mas/pixel and a field of view 
of 0'.'8 x 0'.'8. Further data reduction was done using 
custom-written IDL scripts. For both the science tar- 
get and the telluric calibration stars, we first calculated 
the centroid for each image in the data cubes. Then we 
extracted spectra by using an aperture with a radius of 
10 pixel and a sky annulus between 20 and 25 pixel to 
calculate the flux of the star at each wavelength. The 
spectrum of the telluric standard star contained not only 
the telluric absorption lines, but intrinsic stellar photo- 
spheric lines. To correct for the latter, we used the nor- 
malized solar spectrum available at the SINFONI web- 
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Fig. 1. — Transmission of the atmosphere at Paranal Observatory 
on the night of 30/31 July 2008, in the J, H, and K bands, showing 
lines and bands of telluric absorption. 

pag^E We multiplied it by a blackbody, convolved it to 
the spectral resolution of SINFONI and scaled it. By 
changing the temperature of the blackbody, the width of 
the Gaussian convolution kernel, and the scaling factor, it 
was possible to remove any intrinsic stellar features from 
the telluric correction spectra. The obtained correction 
curves characterize the transmission of the atmosphere as 
a function of wavelength (Fig.[IJ). Since for the first night 
no H and K calibration observations were available, we 
used the average of the curves obtained on the two other 
nights. Finally, the spectra of EX Lup were divided by 
these correction curves, and the results were normalized 
by fitting a third degree polynomial to the continuum. 

Comparison of the spectra taken on the three differ- 
ent nights revealed insignificant differences, thus we av- 
eraged them to increase the signal-to- noise ratio (S:N). 
For most of the following analysis, we used these aver- 
aged spectra, which are plotted in Figs. [21 El and[U S:N 
is 80-160 in the middle of the atmospheric windows and 
far from strong telluric bands (e.g. at 1.23-1.27/im, or 
at 2.21-2.28 /im), 40-80 where there are strong telluric 
absorption (e.g. around 2.01 /im), and 20-40 at the edges 
of the atmospheric windows (e.g. above 1.34 /jm in the J 
band or above 2.40 ^m in the K band). 

2.2. Near-infrared polarimetry 

We have acquired AO-assisted near-infrared polarimct- 
ric observations of EX Lup in visitor mode using the 
NACO instrument mounted on the UT4 telescope of the 
VLT on 10/11 April 2008, as part of the project 381.C- 

0241 (PI: A. Kospal). NACO consists of the Nasmyth 
Adaptiv e Optics System and the CONICA near- i nfrared 
camera (jLenzen et al.l 119981: IRousset et all I2003D . The 
observations were made using the di fferential polarimet - 
ric imaging technique (DPI; see e.g. IKuhn etaTI I2001D . 
The basic idea of the DPI is to take the difference of two 
orthogonally polarized, simultaneously acquired images 
of the same object in order to remove all non-polarized 
light. As the non-polarized light mainly comes from the 
central star, after subtraction only the polarized light, 
such as the scattered light from the circumstellar ma- 
terial, remains. We obtained polarimctric images with 
NACO through the H filter, using a Wollaston prism 

5 These spectra were created from data made available by the 
NSO/Kitt Peak Observatory, produced by NSF/NOAO. 
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with a 2" Wollaston mask to exclude overlapping beams 
of orthogonal polarization, using the 27 mas pixel -1 scale 
camera. EXLup was observed at four different rotator 
angles of 0°, 45°, 90° and 135°. At each angle, a three- 
point dithering was applied. The polarimetric calibra- 
tors were HD 94851 and HD 64299. The data reduction 
was done in IDL using prev iously deve l oped software 
tools present e d in d etails in lApai et all ()2004[ ) and in 
iKospal et al.l (|2008bT ). Calibration measurements of the 
unpolarized standards HD 94851 and HD 64299 indicate 
that the instrumental polarization is less than 3%. 

3. RESULTS AND ANALYSES 

Line inventory. — The near-infrared spectrum of EX Lup 
displays many emission lines. No obvious absorption 
features are visible. We identified the emission fea- 
tures using wavelengths from the NIST Atomic Spectra 
Database^ and fro m the Infrared Spec tral Atlases of the 
Sun from NOAO (jWallace et al_ll99fih . There is no sys- 
tematic difference between the observed and tabulated 
wavelengths; the difference for the individual lines being 
in average less than 2 pixels or 70kms~ 1 . 

The most conspicuous features in the spectrum are the 
lines of atomic hydrogen: the Pa/3 line in the J band 
(the strongest feature in the whole spectrum), the Br 10 
to 20 lines in the H band, and the Br 7 and S lines in 
the K band. We could identify two He I lines: the well- 
known 2.0587 /im feature, and possibly also another line 
at 1.2849 /im (close to the Pa/3 line). Several metallic 
lines are present: lines of neutral Na, Mg, Fe, Si, Ca, Al, 
C, and O are well visible. We could also identify two lines 
of singly ionized Ca, and possibly also the 1.6440 /im for- 
bidden line of singly ionized Fe. Besides the atomic lines, 
molecular bands are also visible: the K band spectrum 
displays pronounced CO ro-vibrational overtone band- 
head emission between 2.29 and 2.40 /im. The species of 
the identified lines are overplotted in Figs. El G2 and|4] 

We analyzed the atomic lines in IRAFu by fitting Gaus- 
sian profiles. We found that the full widths at half max- 
imum (FWHM) of the lines are consistently larger than 
the instrumental FWHM, the latter being 6 A in the mid- 
dle of the J and H bands, and 5 A in the middle of the K 
band, as determined by measuring spectrally unresolved 
OH sky emission lines. Thus, we conclude that the lines 
are spectrally resolved, their measured FWHM are be- 
tween 8 and 12 pixel, and their deconvolved FWHM are 
between 110 and 190 kms -1 . Table [T] shows the tabu- 
lated and observed wavelengths, as well as the equiva- 
lent widths (EWs), of the lines identified in our EXLup 
spectrum. 

Spatial extent. — Since the SINFONI data cubes provide 
spatially resolved information, we checked whether any 
line emission originates from an extended area. We first 
averaged images corresponding to wavelengths in a wave- 
length window of 25 A centered on a certain emission 
line ( "line image" ) . Then we averaged images in a 50 A 
wide window on both sides of the line, at a distance of 

6 http:/ /physics. nist.gov/PhysRefData/ASD/lines_form. html 

7 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



87.5A ( "continuum image"). We created a continuum- 
subtracted line image by simply subtracting the contin- 
uum image from the line image. In order to do PSF- 
subtraction, we supposed that the continuum emission 
is a point source and considered the continuum image 
as a PSF template. Then, we subtracted the PSF tem- 
plate from the continuum-subtracted line image. The 
peak of the PSF template was scaled to the peak of the 
line image. If the line emission is more extended than 
the continuum, the residuals should show this extended 
emission. With this method, we analysed several hy- 
drogen and metallic lines, as well as the CO bandhead, 
but found no extended emission. Thus, we can conclude 
that with our sensitivity, the line emission is not more 
extended than the continuum emission, and that there 
is no extra diffuse line emission apart from the central 
star. The FWHM of the stellar PSF should give a rea- 
sonable upper limit for the size of the emitting region. 
The FWHM was between 4.1 and 6.4 pixels, depending 
on the band and on the efficiency of the AO correction 
for a particular frame. 4 pixels correspond to 50 mas (or 
approximately 8 AU at the distance of 155 pc). Thus, the 
frames with the best AO correction give an upper limit 
of 8AU for the diameter of the near-infrared emitting 
region, both for the continuum and for the H, CO, and 
metallic line emission. A somewhat less stringent upper 
limit can be derived from the FWHM of the broad-band 
H images obta i ned w ith NACO (92 mas or 14 AU). 

lAspin et al.l (|2010t ) observed P Cygni profile for the 
H/3 and the Na D lines in the optical spectrum of 
EXLup obtained in January 2008. The blue-shifted 
absorption i n the se lines appear around — 120 km s -1 . 
I Goto et all (|201lD detected similar blue-shifted absorp- 
tion at — 80kms _1 in the profiles of the v = 1 — fun- 
damental CO lines in August 2008. Outflowing mate- 
rial traveling at a velocity of 80-120 km s -1 would have 
reached an area of 16-24 AU in diameter during the 5 
months that elapsed between the peak of the outburst 
and our SINFONI observations. The fact that we do 
not see emission beyond 8 AU can probably be explained 
by the faint surface brightness of the outflow. This also 
implies that the lines are dominated by emission from in- 
falling or rotating material closer to the star than 8AU. 

Hydrogen lines. — Strengths and profiles of the hydrogen 
lines give information about the physical conditions in 
the emitting medium. In Fig. [5] we plotted the intensi- 
ties of lines in the Brackett series relative to that of the 
Br 7 line. Here we assumed that A y = Omag (see e.g. 
iSipos et al.ll2009l: lAspin et al.ll2010D . and that the shape 
of the SED is flat, thus the absolute line intensities are 
proportional to the line EWs. Then we compared the 
observed intensity ratios to theoretical predictions cor- 
responding to Case B recomb ination (jBaker fe Menzell 
11938- iHummer fe Storevlll987l ). Case B theory is valid if 
the emitting plasma is opaque to Lya photons but op- 
tically thin for higher transition lines. We found that 
our line fluxes are essentially consistent with the Case B 
theory. A model with temperature of 10 000 K and elec- 
tron density of N c =10 7 cm~ 3 best fitted the data points 
(Fig. [5]), however, higher electron densities with higher 
temperatures (N c =10 8 cm -3 , T=12 500K) or lower elec- 
tron densities with lower temperatures (N e =10 6 cm~ 3 , 
T=5000K) also fit within the measurement uncertain- 
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Fig. 2. — Normalized SINFONI J band spectrum of EXLup. Earth symbols indicate telluric absorption. 



ties. The observed Pa/3 / Br 7 ratio of 2.54, however, is 
far from the theoretical Case B value of 6.08. Actually, in 
Case B theory, any combination of temperature between 
1000 and 12 500 K and electron density between 10 2 and 
10 10 cm -3 would give ratios between 4 and 9. This indi- 
cates that the Pa/3 line is probably not optically thin. 

As we mentioned earlier, the hydrogen lines are spec- 
trally resolved. Specifically, the deconvolved FWHM of 
the Pa/3 line is 170±20kms _1 , and of the Br 7 line is 
186 ± 20 km s -1 . Apart from a nearly symmetric broad- 
ened Gaussian peak, high-velocity line wings, extending 
to several hundred km s _1 , can also be seen (Fig. [5]) . Par- 
ticularly interesting is the shape of the Br 7 line wings, 
which is not symmetric: on the blue side, the profile 
agrees with a Gaussian, and starts showing excess from 
— 230kms _1 ; on the red side, the profile agrees with 
a Gaussian up to 340 km s , and shows a pronounced 
bump between 340 and 630 km s -1 . This absence of ex- 
cess emission in the red side may be due to a weak inverse 
P Cygni profile, the signature of mass infall. 

Before the outburst, i n 2001, EXLup di splayed no dis- 
cernible Br 7 emission (Si pos et al.l 120091) . The spectra 
obtained d uring outburst in April 2008 (jJuhasz et al.l 
Isubmittedl ) and in July 2008 (this work) both 
have similar EW (~ — 13 A) and deconvolved FWHM 



(«200kms^ 1 ). Observations by lAspin et all (|2010fl 
show that in February and in May 2008, the EW of the 
Br7 line was slightly larger, —17 A. 

CO bandhead. — In order to model the CO bandhead pro- 
files, we used a simple slab model in which we assume 
that the CO emission comes from a slab of gas in lo- 
cal thermodynamical equilibrium (LTE) with vibrational 
temperature Tco an d CO column density Nco The 
vibrational and the rotational tem perature were set t o 
be equal. We used the equations of IKraus et al.l (fe OOO). 
We to ok the Einstein A coefficients from lChandra et al.l 
(1996), who calculated coefficients up to J=140, and w e 
took the total partition function from IGoorvitchl ([1994). 
We first calculated the absorptio n coefficient p er CO 
molecule, as given by Eqn. 10 in IKraus et al.l ([2000), 
using a Gaussian CO line profile with a width of 2 
km s -1 , as s uggeste d for the turbulent broadening by 
iNaiita et al.l ()1996af ). We only considered temperatures 
between 2000 K and 5000 K. For LTE, we need suffi- 
ciently high temperature and sufficiently high density so 
that the vibrational levels can be collisionally excited. 
This requires temperatures Tco >2000K and densities 
greater than nn > 10 10 cm~ 3 ; but temperatures should 
also be les s than 5000 K, other wise CO molecules would 
dissociate (Scoville et al.|[l980l) . Then, we convolved the 




absorption coefficient as a function of waveiength with 
the velocity profile of a disk in Keplerian rotation (see 
below). We then multiplied the absorption coefficient by 
the column density of CO in order to obtain the optical 
depth. From the optical depth, we calculat ed line inten- 
sities using the transfer equation (Eqn. 9 in lKraus et al.l 

lAspin et al.l (|2010D needed to introduce a velocity pro- 
file of a disk in Keplerian rotation in order to fit their 
observations of the CO bandhead emission. They found 
that the best fit could be obtained using a stellar mass of 
M*=0.6M Q , inner disk radius of r; n =0.08 AU, outer disk 
radius of r ou t=0.13 AU, temperature of Tco=2500K and 
optical depth tqq < 0.1. We calculated the line profile 
emerging from such a disk, assuming a disk inclination 
of i=45° and a disk ra dial brightness p rofile proportional 
to r~ 2 - 5 (same as in iGoto et all 120111) . The resulting 
line profile is plotted in the small inset in Fig. [7] We 
used this profile when calculating the optical depth (see 
above). Finally, we smoothed our model spectrum to the 
instrumental resolution of SINFONI. The result, plotted 
with a gray line in Fig. fits our observations of the 
v=2— >Q, v=3— >1, and v=4— >2 transitions very well. We 
note that a model with an oute r radius of r out =0.4AU, 
as derived by IGoto et al.l (|2011l ) from the profiles of the 
CO fundamental emission lines of EXLup, would fit our 



observations equally well. 

Metallic lines. — M-type stars display numerous Na, Ca, 
K, Fe, Ti, and Si absorption lines in their near-infrared 
spectra , which are often us ed for spectral type classifi- 
cation (|Cushing et al.ll20 05). In quiescence, EXLup also 
displayed atomic absorption lines (SI at 1.199 (im, Mg l 
at 1.577/im, and All at 1.673 fim, see lSipos et al.i r2009). 
The se lines are most pr obably photospheric in origin (see 
also lHerbig et al.ll2001| ). During outburst, no absorption 
lines are visible any more, despite the fact that, at the 
time of our SINFONI observation, 77% of the J and ri- 
band flux, and 60% of the K-band flux was photospheric. 
These values were obtained by supposing that the quies- 
cent JHK photometry represent the stellar photosphere, 
and that during our observations, the V— J, V— H, and 
V— K colors were the same during ou r observations as in 
April 2008 (|Juhasz et alJ lsubmittcd). In outburst, the 
above mentioned Mgl line is in emission, the lines that 
were in emission already in quiescence (the Pa (3 line and 
an O I line at 1.129 fim.) became much stronger, and many 
more atomic emission lines appeared. 

Most of the metallic lines visible in the SINFONI spec- 
trum originate from neutral atoms whose ionization po- 
tential is low (in the range of 5.1— 8.2 eV, cf. 13.6 eV for 
hydrogen), the only exception is oxygen (with an ioniza- 
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Fig. 4. — Normalized SINFONI K-band spectrum of EXLup. The earth symbol indicates telluric absorption. 
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tion potential of 13.6 eV). This means that these atoms 
are located in an area where hydrogen is mostly neu- 
tral, probably shielded from most of high energy radia- 
tion from the accretion shock region by the circumstel- 



lar disk itself. However, the atoms cannot be located 
too deep in the disk, since many lines and line ratios 
indicate fluorescent emission pumped by UV photons. 
The Nal doublet at 2.206 and 2.209 /im may be pumped 
by 3303 A continuum photons. If this is the case, then 
the Nal line at 1.141 /im should also be present with a 
strength roughly twice of those of the lines at 2.206 and 
2.209 /im, while collisional excitation would cause a line 
ratio above 10 ([McGregor et al.lll988l ). For EXLup this 
ratio is wl.3, supporting the fluorescent origin. Another 
fluorescent line is that of OI at 1.129 /i m, which may be 
pump ed by Ly /3 or continuum photons (|McGregor et alj 
119841 and references therein). Continuum fluorescence 
predicts an OI 1.316 /im / 1.129 /im line ratio above 
unity. The 1.316 /xm line is practically invisible in our 
SINFONI spectrum, which favors the Ly /3 fluorescence. 
If this is true, another O I line at 8446 A should also be 
present. Indeed, this line is well visible in our FEROS 
spectra of EX Lup obtained be tween April and June 2008 
(jSicilia-Aguilar et al.llin prep.l ). Interestingly, both the 
8446 A and the 1.129 /im lines are visible even in quies- 
cence, although they are much weaker than in outburst, 
indicating that a sufficient amount of Ly /3 radiation was 
present even in quiescence. The presence of Ly /3-pumped 
lines of neutral oxygen indicates the presence of dense 
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Fig. 6. — Solid black line: profile of the Pa and Br 7 lines; 
dashed gray line: fitted Gaussian profile. 

warm regions embedded in hot ionized plasma. The Mg I 
1.503 /mi line may also be formed by UV fluorescence. 
This requires that the Mgl 1. 574/1. 575/1. 577/im transi- 
tions are also present. This is true for EXLup. However, 
collisional excitation may also play a role here, since the 
Mgl 1.183 /im line is also visible. 

Shock signatures. — We searched the spectrum of EX Lup 
for shock-excited lines. We found that the H2 line at 
2.1218 /zm is not visible. The [Fell] line at 1.6440 jum 
may be present, but its unambiguous identification is 
difficult due to its closeness to the strong Br 12 line. We 
checked the position-velocity diagram and the spectro- 
astrometric signal (see below) in order to identify the 
location of the 1.6440 /jm emission. We found neither 
extended emission nor positional offset compared to the 
continuum position. Thus, supposing that the [Fell] line 
is indeed present, at the time of our observations the 
shock was still very close to the central source. This 
apparent lack of strong forbidden lines may be an indi- 
cation that the emission line region has a high densit y 
where forbidden lines are quenched (|Nisini eTaD[2005). 
Further monitoring of EX Lup may reveal shocks as they 
reach less dense regions farther from the central star, ex- 
panding into the medium around the disk. 

Spectro-astrometry. — The SINFONI data cubes can be 
used for spectro-astrometry, i.e. measuring the position 
of the source as a function of wavelength. If there is 
extended emitting material moving at different veloci- 
ties, the source position measured at these velocities will 
deviate from the position measured at continuum wave- 
lengths. In case the emission comes from a rotating disk, 
at certain wavelengths, there will be a positive offset in 
the centroid, while at other wavelengths, there will be a 
negative offset. Thus, depending on whether we see the 
approaching or the receding part of the disk, the mea- 
sured spectro-astrometric signal as a function of velocity 
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wi ll follow a typical sinuso idal pattern (see e.g. Fig. 2 
in iPontoppidan et al.l 12008). Infall or outflow would, in 
most cases, cause a positive only or neg ative only signal 
(see the modeling bv lEisner et alJl2010j ). For certain ge- 
ometries, a bipolar outflow could cause both positive and 
negative offsets for different velocities, but the centroid 
displacement w ould be more linear (see e.g. the example 
of V536Aql in iWhelan et al.l 12004ft . Binaries may also 
cause a non-zero spectro-astrometric signal, but again, a 
positive only or negati ve only signal (|Takami et al.l[2003t 
iBrannigan et al.l 120061 ). 

We analyzed several emission lines in the spectrum of 
EXLup with spectro-astrometric technique. Fig.|5]shows 
some examples. We found that only the hydrogen lines 
indicate any spectro-astrometric signal. Metallic lines 
show a flat signal, which can be either due to gas close to 
the star, or slow-moving gas at larger distances. The fact 
that most of the metallic atoms are neutral supports the 
second possibility. The spectro-astrometric signals plot- 
ted in Fig.[8]for the Pa j3, Br 7, and Br 11 lines are all indi- 
cating rotating material. Since SINFONI contains no slit, 
our measurements do not suffer from the artifacts usual 
for slit spectra, and the spectro-astrometric signal can be 
calculated for any angle across the source. The results 
can be then used to calculate the position angle of the 
rotation axis, i.e. the angle where the spectro-astrometric 
signal is constant 0. In our case, this yields P.A. = 
80 ±10°, thus, the rotation axis is almost exactly east- 
west oriented. Remarkably, the spectro-astrometric sig- 
nals for the hydrogen lines are all symmetric within the 
measurement uncertainties, indicating that the bright- 
ness and velocity distribution of hydrogen gas around 
the central star is azimuthally symmetric. 

Polarimetry. — Broad-band polarfmetric images through 
the H filter were obtained for EXLup using a shorter 
(0.35 s) and a longer (2 s) exposure time. The short ex- 
posure images show that the distribution of the total in- 
tensity follows a Gaussian profile. The profile of the long 
exposure images indicates that the inner region (within 
a radius of 3 pixels) of the source is saturated, but the 
rest of the profile is the same as for the short exposure 
images, suggesting no extended emission. The Stokes Q 
and U images, as well as the polarized intensity images 
show no polarized light whatsoever. This limits the H- 
band degree of polarization for EX Lup below the instru- 
mental polarization (not more than 1-2%). This result 
is not surprising if we consider that only a small fraction 
of the observed H-band flux is expected to be scattered 
light (thus polarized) ; most of it is photospheric emission 
from the central star, or thermal emission from the disk. 

4. DISCUSSION 
4.1. The origin of the CO emission 

Our model indicates that the CO bandhead emis- 
sion is consistent with emission from CO gas which has 
a temperature of 2500 K and which is located in the 
disk between 0.08 and 0.13 AU from the central star. 
Uuhasz et al.l ([submitted! ) found that in outburst there is 
a dust-free inner hole in the system within 0.3 AU. Thus, 
the CO gas we observe here should mostly be located in 
t his dust-free region . 

iGoto et ahl (|2011f ) also derived the properties and dis- 
tribution of the CO gas using the fundamental vibra- 
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Fig. 7. — CO overtone bandhead emission in the spectrum of EXLup. The gray line indicates a model of 2500 K CO gas, convolved 
with a line profile characteristic of a Keplerian disk around a 0.6 Mq star with inner and outer radii of 0.08 and 0.13 AU, respectively, and 
inclination of 45° (see small inset). 



tional lines in the 4.6-5 /im wavelength range. They 
found that the fundamental lines are the superposition of 
a narrow (FWHM=50 km s" 1 ) and a broad (FWZI=150 
km s _1 ) component. The broad component comes from 
gas orbiting the central star at 0.04-0.4 AU, and has a vi- 
brational temperature between 1800 and 3000 K. These 
parameters are very similar to those we obtained from 
the overtone bandhead features. Thus the gas emitting 
the overtone and the fundamental lines are co-located 
and are very likely physically the same material. 

The CO bandhead feature in EXLup displays signifi- 
cant t emporal variability: it was in absorp tion in quies- 
cence (|Sipos et al.H2009t IBerbig et al.ll200lD . very strong 
emission in Febru ary 2008, weak emission in May 2008 
(|Aspin et al.ll2010[). and also weak em ission in July 2008 
(this work). lLorenzetti et al.l (|2009f) observed the vari- 
ability of the CO bandhead feature in several other EX- 
ors. They claim that CO absorption is associated with 
quiescent periods, when the accretion rate is low and the 
stellar photosphere (of an M-type dwarf) is visible. Ac- 
cording to them, CO emission is associated with more 
active periods, when the accretion rate is higher, UV ra- 
diation is higher as well, thus the CO gas in the inner disk 
is heated and is producing CO bandhead emission. The 
general behavior of the CO bandhead feature in EX Lup 
is consistent with this picture. 

There are several scenarios described in the literature 
to explain the origin of CO bandhead emission in YSOs. 
Based on observations of low-mass YSOs indicating sig- 
nificant CO bandhead e mission variab i lity o n timescales 
as short as a few days, IBiscava et al.l (|1997l ) considered 
the disk, stellar/ disk wind, and funnel flows as the origin 
of CO emission. Although our observations are consis- 
tent with the inner disk origin (Fig. [7]), we cannot exclude 
the funnel flow origin either: the figure in [Martini (|1997D 
suggests that funnel flows can produce CO bandhead pro - 
files very similar to those we observed. IGoto et"aLl (|2011[ ) 
argued against the funnel flow origin based on the lack 
of systematic redward absorption in the fundamental CO 
lines. However, one can imagine a situation when - due 
to the inclination of the system - the funnel flows do not 
cross the line-of-sight toward the hot spot on the stellar 
surface, thus, absorption is not expected to be seen. We 
note that CO can survive without dissociation only in the 
outer parts of the funnel flows, close to the inner edge of 



the disk, thus, there is no essential physical difference 
between a funnel flow and an inner disk origin. 

In order to test whether the heating of the CO gas 
is related to accretion, at different epochs we correlated 
the strength of the CO bandhead feature with the visual 
brightness of the star (which is supposed to be accretion 
luminosity) and with the Br 7 line flux (another usual ac- 
cretion rate tracer). Using t he spline-smoothed versio n 
of the light curve presented in lJuhasz et al.l (jsubmittedl ) . 
we estimated a visual brightness of 9.5 mag in February, 
10.2 mag in May, and 10.2 mag in July 2008. Assuming 
that the V— K color of the source was similar at these 
epochs to the value me asured in April 2008 (2.05 mag, 
Uuhasz et alJIsubmittedl ). the Br 7 line fluxes decreased 
by a factor of 1.9 between February and May, and a fac- 
tor of 2.3 between February and July. We calculated 
the EWs of the CO bandhead by integrating the spectra 
between 2.2925 and 2.2975 /im, and obtained —21 A for 
February, -4 A for May, and -3.3 A for July 2008. As- 
suming again a constant V— K color, the flux of the CO 
bandhead decreased by a factor of 10 between February 
and May, and a factor of 12 between February and July. 
Thus, the strength of the CO bandhead feature, the op- 
tical brightness of the star, and the strength of the Br 7 
feature seems to be changing in the same direction, al- 
though not exactly with the same rate. Consequently, 
our results indicate that both heating by irradiation and 
heating by accretion may have a role in the excitation of 
the CO ban dhead. We note t hat the broad component 
observed bv lGoto et al.l (|2011l ) in the fundamental emis- 
sion lines was also fading with time, but due to different 
observing epochs, we cannot directly compare the results 
for time evolution of the fundamental and the overtone 
lines. 

4.2. The origin of metallic lines 

The lack of spectro-astrometric signal for the metal- 
lic lines and the fact that most of the metallic lines are 
from neutral atoms indicate that these lines are originat- 
ing from regions farther from the central star than the 
hydrogen lines. Some of the lines are fluorescently ex- 
cited by UV photons, indicating that at least part of the 
atoms are subject to direct stellar radiation, which con- 
dition is probably fulfilled at least on the surface of the 
inner gas disk. Thus, the metallic atoms may be partially 
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Fig. 8. — Top panels: Profile of the Pa/3, Br7, Br 11, and the 
1.4882 fim Mgl line (the flux is normalized to the continuum). Bot- 
tom panels: Displacement of the of the image centroid measured on 
the line+continuum images with respect to the image centroid at 
continuum wavelengths. The model overplotted with a gray line on 
the Br 7 spectro-astrometric signal is a Keplerian disk model. For 
model parameters, and a more detailed discussion of the spectro- 
astrometric data, see text. 



co-located with the CO gas. Interestingly, the 2.206 and 
2.209 /im Na doublet shows a similar time evolutio n to 
that of the CO bandhead. lLorenzetti et al.l (|2009t ) no- 
ticed that it follows the same behavior as the CO band- 
head in several other EXors: either both are in emission 
or both are in absorption. They interpreted this result 
as a common origin: both the Na and the CO feature is 



originated in the stellar photosphere when in absorption, 
an in the inner disk when in emission. Our results for 
EXLup is consistent with this s cenario: both fea tures 
were in absorption in quiesce nce (ISipos et al.ll2009[) , and 
were in emission in outburst (|Aspin et al.H2010i and this 
work) . 

4.3. The origin of the hydrogen emission 

Hydrogen emission in YSOs is generally believed to be 
origi nating in disks, funnel flows or stella r /disk winds 
(e.g. iNaiita et al1l!996bt iKraus et all 120101) . In the fol- 
lowing, we check whether any of these scenarios are con- 
sistent with our results on EXLup. The excitation di- 
agram in Fig. [5] indicates that the hydrogen is hotter 
than the CO gas, suggesting the hydr ogen gas is closer to 
the st ar th an 0.04-0.08 AU. Fo llowing iPontoppidan et al.l 
(120081) and lEisner et al.l l|2010l) . we modeled the Br 7 line 
using this scenario and adopting a simple disk model. We 
supposed that all continuum arises from (and only from) 
the central star, and all the Br 7 emission originates from 
(and only from) the disk. We assumed that the velocity 
field in the disk is Keplerian around a 0.6 MQ-mass cen- 
tral star, while the brightness profile is a power law. With 
these assumptions, we calculated the line-of-sight com- 
ponent of the velocity, and made images corresponding 
to each SINFONI velocity channel. We added a central 
point source with a wavelength-independent brightness, 
and blurred the image by convolving it with a Gaus- 
sian PSF. We measured the spectro-astrometric signal 
on these synthetic images the same way as we did for 
the real observations. We used a disk extending from 
the stellar surface (0.007 AU) to 0.04 AU, inclination of 
45°, and a radial brightness profile proportional to r~ 2 - 5 . 
The resulting line profile and spectro-astrometric signal 
are overplotted with gray lines in Fig. [5] The results 
indicate that emission from such a Keplerian disk can- 
not significantly contribute to the observed Br 7 line, be- 
cause the emerging line profile would be double-peaked 
and the model spectro-astrometric signal indicates that 
high- velocity gas is closer to the star than observed. 

We calculated different disk models by changing the 
disk inclination, disk inner and outer radii, and power 
law exponent of the disk brightness profile. We could not 
find a model that fits the observed spectro-astrometric 
signal. This analysis seems to exclude the equatorial 
boundary layer as the main path of accretion onto the 
star, a model often invoked for FU Orionis-type out- 
bursts (jHartmann fc Kenvonl [19961 ) . especially since the 
innermost parts of a boundary layer is sub-Keplerian. 
Although the boundary layer, or strong turbulence 
could explain the single-peaked line profile, they can- 
not reproduce the high- velocity high-amplitude spectro- 
astrometric signal. 

Our results suggest that hydrogen gas is orbiting the 
central star but its velocity profile is not Keplerian, in 
a sense that there is high-velocity gas farther from the 
star. This points to a funnel flow or disk wind origin, 
where the material is launched from the inner disk along 
the magnetic field lines with high velocities. This may 
also help to explain why the hydrogen gas is hotter than 
the CO. The wind scenario is supported by the P Cygni 
profile s observed in several optical lines by lAspin et al.l 
(2010). Supposing that the disk wind scenario is true 
for EXLup, it should be optically thin, otherwise the 
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excitation diagram in F ig. [5] would show natter line ratios 
(|Lorenzetti et al.ll2009t) . 

5. CONCLUSIONS 

In this paper we performed a near-infrared spectro- 
scopic and spectro-astrometric study of EXLup during 
its most extreme outburst in 2008. Our main conclusions 
are the following: 

• The CO bandhead emission comes from an in- 
ner, dust-free region. This area is id entical with 
the br oad component area identified bv lGoto et al.l 
(|2011l ). from where the CO fundamental emission 
comes from. The CO gas is rotating around the 
central star, either at the inner edge of the disk or 
in the outer parts of funnel flows. 

• EXLup also displays numerous neutral atomic 
lines, which probably come from an area similar 
to that of the CO gas. At least part of this area is 
subject to direct UV photons, as evidenced by the 
presence of fluorescent emission lines. 

• A very conspicuous feature in the spectrum of 
EXLup is the Brackett series, whose excitation di- 
agram indicates that the hydrogen emission is op- 
tically thin. Spectro-astrometric analysis of these 
lines suggests that the hydrogen emission is proba- 
bly not coming from an equatorial boundary layer; 
a funnel flow or disk wind origin is more likely. 

Based on our findings, we can attempt to reconstruct 
the geometrical and kinematic structure of the circum- 
stellar material within the inner few tenths of AU in 
the EXLup system during the outbu rst. The dust disk 
has an inner radius o f 0.2-0.3AU (jSipos et al.l 120091 : 
Uuhasz et alJIsubmittedl ). Within this area, there is an 
optically thin gas disk, whose temperature is a few thou- 
sand K, as indicated by the CO and neutral metal emis- 
sion lines. This area certainly also contains hydrogen 
gas, but this component is not visible in the spectro- 
astrometric signal due to its low velocity. We see ev- 
idence for high temperature (psIOOOOK), high velocity 
(wlOOkms" 1 ) hydrogen gas in the system, which is not 
located in the equatorial plane. Part of this hydrogen 
probably falls onto the stellar surface along magnetic 
funnel flows. Inde ed, based on X-ray and UV data, 
iGrosso et al.l (|2010T ) reported on the presence of accre- 
tion shocks and accretion hot spots on the stellar sur- 
face of EX Lup. Some of the UV photons thus generated 
must reach the disk surface and produce the fluorescent 
emission we observe in certain sodium, magnesium and 
oxygen lines. The UV radiation may also trigger chem- 
ical changes in the disk. Part of the hydrogen gas does 
not fall onto the stellar surface but leaves the system in 
the form of a hot wind, as evidenced by the P Cygni 
profile of the optical hydrogen lines. Considering all 



these arguments, the emerging picture is broadly consis- 
tent with that of the standard magnetospheric accretion 
model us ually assumed for n ormally accreting T Tauri 
stars fe.g. IBouvier et £tl.|[2007h . 

Several models are described in the literature to ex- 
plain the increased accretion in young eruptive stars. 
Our results may place constraints on the applicability of 
these models for the case of EX Lup. Our results do not 
indicate the presence of a fully ionized ring of material 
during the eruption of EX Lup, whi ch seems to contradict 
the th ermal instability model of lHartmann fe Kenyonl 
(|1996[ ) and iBell fe Linl (|1994h . No stellar or plane- 
tary companions to EXLup have been found so far 
(jSipos et al.ll2009l and references therein), neither do our 
spectro-astrometric analysis suggests the presence of one. 
This makes outburst mod els involving perturbatio n by a 
close co mpanion unlikely dBonnell fe Bastier] [l992'l. The 
model of iVorobvov fe Basul (|2010f ) involves gravitational 
instability and fragmentation in the outer disk, and the 
infall of these fragments onto the star. The modest disk 
mass, and the fact that our spectro-astrometric observa- 
tions indicate an azimuthally symmetric mass distribu- 
tion in the inner disk reasons against this model. 

Recently. iD'Angelo fe Spruit] (|2010l ) proposed that ac- 
cretion onto a strongly magnetic protostar is inherently 
episodic if the disk is truncated close to the corotation 
radius. In their model, the magnetic field initially trun- 
cates the disk outside the corotation radius, thus accre- 
tion onto the star is inhibited. As gas in the inner regions 
of the disk piles up, material is pushed inside the coro- 
tation radius, and the accumulated material is accreted 
onto the star until the reservoir is depleted, and the inner 
radius of the disk is again outside the corotation radius. 
ISipos et al.l (|2009| 1 reported a v sin i of 4.4 ±2 km s 1 for 
EXLup, which translates into a rotation period of about 
13 days (using an inclination of 45° and stellar radius of 
1.6 Rq), and a corotation radius of about 0.3AU (using 
a stellar mass of 0.6 M Q ). This value is close to the ra- 
dius of the dust-free zone. Our results show that during 
outburst, gas is present in this area, potentially chan- 
neled along magne t ic fun nel flows. Thus, the model of 
ID'Angelo fc Spruit] EttWh may be applicable for EX Lup. 
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TABLE 1 

Spectral features in the near-infrared spectrum of EX Lup. 



Species 


\ab [^ m ] 


^obs \p>m] 


EW [A] 


Comment 


Fe I 


1.1123 


1.1131 


—0.73 




Na I 


1.1200 


1.1200 


—0.44 




A 1 T 
Al 1 


1.1258 


1.1261 


—0.52 




O I 


1.1290 


1.1293 


—7.97 




IN a 1 


1.1385 


1.1386 


— 1.19 




Na I 


1.1407 


1.1409 


— 1.60 




Fe I 


1.1425 


1.1428 


—0.45 




T7V, T 

re 1 


1.1442 


1.1445 


—0.61 




be 1 


1.1597 


1.1599 


—0.93 




Fe I 


1.1610 


1.1613 


— 1.69 




Fe I 


1.1641 


1.1644 


— 1.23 




t 

Kj 1 


1.1662 


1.1665 


—0.38 




Fe I 


1.1693 


1.1696 


— 1.30 




t 

Kj 1 


1.1757 


1.1758 


— 1.57 




Fe I 


1.1786 


1.1789 


— 1.16 




Mg I 


1.1831 


1.1834 


—2.02 




/~1 _ TT 

Ua 11 


1.1842 


1.1845 


—2.17 




re 1 


1.1886 


1.1889 


—3.30 




C I 


1.1899 


1.1901 


—0.30 




Ca II 


1.1953 


1.1955 


— 1.58 




T 

re i 


i . ±y t d 


1 1 Q7Q 

i . ±y t y 


— z.oo 




Q; t 

Ol 1 


i.iye f 


1 i onn 

i.iyyu 


1 A Q 

— 1.4o 




Qi T 
Ol 1 


1 i one; 

i.iyyo 


1 1 007 

i.iyy ( 


ill 
— 1.11 




Qi T 
ol 1 


l.ZUoO 


1 OOQ7 

l.ZUo I 


O 1 o 
— Z. 1Z 




Mg 1 


i onc7 
l.ZUo I 


i.zuoy 


n 71 
— U. / 1 




til T 
ol 1 


1 01 D7 
1.Z1U I 


1 01 no 
i .ziuy 


n n/i 
— U.y4 




O; T 
Ol 1 


l.zz / 4 


1.2276 


n qo 
— U.oz 




H I 


1.2821 


1.2824 


—29.34 


Pap 


TT „ T 

HC 1 


i oo/i n 
l.Zo4y 


1 OQ£/1 
l.Zo04 


n on 
— U.ZU 




A 1 T 
Al 1 


L.olZ ( 


1 q 1 on 

i.oizy 


n cri 
— U.01 




A 1 1 
Al 1 


l.olD4 


1 Q 1 X."7 
l.olO ( 


— U.oo 




TjV. T 
re 1 


l.ozyz 


i.ozyo 


n a 1 
— U.41 




Fe I 


1.4547 


1.4550 


— 1.82 




Mg I 


1.4882 


1.4884 


—2.75 




Mg I 


1.5029 


1.5032 


—4.49 




Mg I 


1.5044 


1.5047 


—3.39 




Mg I 


1.5052 


1.5055 


— 1.76 




T T T 

H 1 


1.5195 


1.5200 


—0.35 


Br 20 


TT T 

H 1 


1.5265 


1.5267 


—0.36 


Br 19 


Fe I 


1.5299 


1.5300 


—0.24 




H I 


1.5346 


1.5347 


—0.82 


Br 18 


H I 


1.5443 


1.5447 


—0.63 


Br 17 


TJ T 

H 1 


1.5561 


1.5562 


— 1.11 


Br 16 


T7V, T T-J T 

re 1, rl 1 


l.ooy7,1.57Uo 


1.5708 


1 AT 

— 1.47 


T7V-, T 1 D. 1 C T 1 _„ J 

re 1 + r>r lo blend 


Mg I 


1.5745 


1.5747 


—0.38 




T\ T _ T 

Mg 1 


1.5753 


1.5757 


—0.99 




Mp- T 


1.5770 


1.5774 


— 1.51 




H 1, iVlg 1, Ol 1 


i.oooo, i.ooyi, i.ooyo 


i.Ooyz 


Q OA 
— 0.Z4 


rsr 14 -+- ivig 1 oi 1 blend 


T 

Ol 1 


i.oyo4 


1 c;nfi7 

i.oyo ( 


n /17 
— U.4 / 


1 

Ol 1 


l.uUyy 


1 ai no 
1.01UZ 


n kq 
— U.OO 




tt t 
rl 1 


1 C11/| 

1 .01 14 


1.011 i 


O 1 7 
— Z. 1 ( 


D„ 1 Q 
DI lO 


CS T TJU T 

ol 1, re 1 


l.OooO, l.Ooof 


l.OoyU 


n 7Q 
— U. / o 


blend 


TT T 

rl 1 


1.041Z 


1.0410 


O 70 

— Z. / z 




fnv^ tti 
[Fc Ilj 


1.6440 


1.6436 


n *? a 
— U. (4 




H I 


1.6811 


1.6815 


-3.56 


Br 11 


C I 


1.6895 


1.6900 


-0.39 




Mg I 


1.7113 


1.7117 


-0.72 




H I 


1.7367 


1.7369 


-4.36 


Br 10 


H I 


1.9451 


1.9455 


-8.87 


Br 8 


Cal 


1.9511 


1.9515 


-1.20 




Cal 


1.9782 


1.9786 


-2.43 




Cal 


1.9868 


1.9871 


-0.98 




He I 


2.0587 


2.0589 


-0.90 




H I 


2.1661 


2.1665 


-11.54 


Br 7 


Nal 


2.2062 


2.2066 


-1.13 




Na I 


2.2090 


2.2093 


-1.32 





